arXiv:1505.01272vl [cond-mat.quant-gas] 6 May 2015 


Self-Accelerating Matter Waves 

C. Yuce 

Department of Physics, Anadolu University, Turkey □ 

(Dated: May 7, 2015) 

The free particle Schrodinger equation admits a non-trivial self-accelerating Airy wave packet 
solution. Recently, the Airy beams that freely accelerate in space was experimentally realized in 
photonics community. Here we present self-accelerating waves for the Bose-Einstein condensate in 
a time dependent harmonic oscillator potential. We show that parity and time reversal symmetries 
for self accelerating waves are spontaneously broken. 


I. INTRODUCTION 


In 1979, Berry and Balazs theoretically showed that 
the Schrodinger equation describing a free particle ad¬ 
mits a non-trivial Airy wave packet solution [l|. This 
free particle wave packet is unique in the sense that it 
accelerates. Furthermore, the Airy wave packet doesn’t 
spread out as it accelerates. The Airy wave packet is also 
called self-accelerating wave packet since it accelerates 
in the absence of an external potential. The accelerat¬ 
ing behavior is not consistent with the Ehrenfest theo¬ 
rem, which describes the motion of the center of mass of 
the wave packet. The reason of this inconsistency is the 
non-integrability of the Airy function. The Berry and 
Balazs’ paper initiated extensive theoretical investiga¬ 
tion on such waves [2^ • The self-accelerating Airy wave 
packet was also experimentally realized within the con¬ 
text of optics three decades after its theoretical prediction 
[in. Since then, self-accelerating wa ve p ackets have 
stimulated growing research interest |15l422j| . Soon after 
the first experimental realization of the self-accelerating 
Airy beams in optics, the same types of beams using free 
electrons instead of light was generated [l^. Using a 
nanoscale hologram technique, the Airy wave packet for 
electrons is obtained by diffraction of electrons. The self- 
accelerating and non-spreading solution is not restricted 
to the linear Schrodinger equation and nonlinear gener¬ 
alization was considered by some authors [l^, [HI . 

The ultracold systems are used to study the effects in 
quantum and statistical physics both theoretically and 
experimentally. Experimental achievements and theoret¬ 
ical studies of Bose Einstein condensates (BEG) of weakly 
interacting atoms have stimulated intensive interest in 
the field of atomic matter waves [H, [2^. Recently, self 
accelerating matter waves are considered in the absence 
of external potential [2^. In experiments of BEG, har¬ 
monic potential is mostly used to trap atoms. Therefore 
self accelerating matter waves in the presence of external 
harmonic potential is worth studying. In this paper, we 
find and examine self accelerating wave packets for ultra¬ 
cold atoms in a time dependent harmonic trap. We derive 
a formula for self-acceleration of such waves and show 
that the self-acceleration depends on the initial form of 
the wave packet. 


II. SELF ACCELERATING AIRY WAVES 


In the present study, we consider ultracold atoms in 
a time dependent trap, which can be realizable experi¬ 
mentally [^. The physics of ultracold gases in a time 
dependent trap was theoretically investigated more than 
a decade ago jH, [2^. Another physical freedom on the 
study of ultracold gases is the temporal tunability of non¬ 
linear interaction. Some experiments on ultracold experi¬ 
ments have demonstrated that tuning of the nonlinear in¬ 
teraction strength can be achieved by applying an exter¬ 
nal magnetic field, known as the Eeshbach resonance [30| . 
With these experimental degrees of freedom, it is possi¬ 
ble to study self-accelerating waves for ultracold atoms. 
The physics of ultracold atoms is well described by the 
Gross Pitaevskii (GP) equation [3l|. The 1-D GP equa¬ 
tion for a BEG with time dependent nonlinear interaction 
strength in a time dependent harmonic trap reads 






dt I 2m dx^ 




iIj ( 1 ) 


where oj{t) and g(t) are the time dependent angular fre¬ 
quency and nonlinear interaction strength, respectively. 
We assume that the nonlinear interaction strength is pos¬ 
itive for all time. In the non-interacting limit, the system 
obeys the Ehrenfest’s theorem provided that the wave 
function is square integrable. According to the Ehren¬ 
fest’s theorem, the acceleration of wave packet can be 

^ X ^ 

found using the equation- — -h < x >= 0. If 

dP 

either the wave packet is initially displaced from the ori¬ 
gin, i.e., < x(0) 0 or the initial velocity of the wave 

packet is different from zero, i.e. < x(0) >= 0, then the 
wave packet accelerates. Although the additional pres¬ 
ence of the nonlinear interaction changes the density pro¬ 
file, it has nothing to do with the acceleration of the wave 
packet. The Ehrenfest’s theorem shows that the expecta¬ 
tion value obeys the classical dynamical laws. Therefore 
the acceleration of the wave packet matches the accelera¬ 
tion of the classical particle moving in one dimension un¬ 
der the influence of harmonic potential. This statement 
is true only when the wave packet is square integrable. 
Below we will find non-integrable wave packet solution 
of the GP equation. In other words, we are looking for a 
wave packet that doesn’t obey the Ehrenfest’s theorem. 
To get self-accelerating wave packet solution, let us first 
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rewrite the GP equation in an accelerating frame, 


, X-Xc{t) 

X = -- 

L{t) 


( 2 ) 


where the time dependent function Xc{t) describes trans¬ 
lation and L{t) is a time dependent dimensionless scale 
factor to be determined later. More precisely, we will see 
that the width of the wave packet changes according to 
Lit). Initially, we take i(0) = 1. Under this coordinate 
transformation, the time derivative operator transforms 
as dt ^ dt — {L x' + Xc)/L dx', where dot denotes time 
derivation. 

In the accelerating frame, we will seek the solution of the 
form 




( 3 ) 


where the position dependent phase reads 
A{x',t) = ^ (ax' + ^x'^ + Sj and the time de¬ 


pendent functions are given by a{t) = Lxc, /3{t) = LL 
1 

and Sit) = -Xc^ -Substitute these transfor- 

2 2 

mations into the GP equation and assume that the 
following equations are satisfied 


L + uj^it) L = 0, (4) 

Xc+Uj'^it)Xc = (5) 


where Og > 0 is an arbitrary constant in units of acceler¬ 
ation that can be experimentally manipulated as we shall 
see below. 

Suppose the nonlinear interaction changes according to 
git) = go/Lit), where go > 0 is a constant. Then the 
time dependent GP equation is transformed to the time 
independent second Painleve equation [s^ 


2m^ao 


dx'^ 


x' + GolV’P W” = 0. 


( 6 ) 


where Go = 2m/‘h^gQ. 

One can find a non-integrable stationary solution of 
this equation. By transforming backwards, the wave 
packet in the lab. frame can be obtained. However, 
such a non-integrable solution can’t be used to find the 
acceleration of the wave packet, < x >. Fortunately, 
the equations glSl) describe time dependent width 
and acceleration of the original wave packet. Before 
getting stationary solutions of the second Painleve 
equation, let us discuss these equations in more detail. 
The equation for L(<) implies that spreading of the 
wave-packet is uniquely determined by the external time 
dependent harmonic potential. The wave packet remains 
non-dispersive in the absence of external potential, 
w = 0. Depending on the angular frequency a;(t), one 
can also obtain breathing or expanding/contracting 
wave packet solution. The equation m implies that 


the motion of the wave packet depends not only on 
the external harmonic potential but also on the scale 
factor Lit) and the constant og. The constant og plays 
a vital role for self accelerating waves since wave packets 
follow the trajectory of the classical motion in the limit 
ag = 0. The self-acceleration depends also on Lit), 
which is determined by the external potential. The 
presence of external potential expands/contracts the 
self accelerating wave and the expansion/contraction 
contributes acceleration. It is interesting to note that 
an initial phase of the wave function plays a role on 
the self-acceleration. As can be seen from the equation 
o, different initial conditions T(0) lead to different 
scale functions L(t) and this contributes self acceleration 
according to ([S]). In other words, imprinting the phase 
on the initial wave function, i.e., /3(t = 0) 0, changes 

the self-acceleration of the wave packet. 

Let us now discuss the solutions of the equations gia 
for some certain cases. The self accelerating solution 
obtained by Berry and Balasz is recovered when to it) = 0 
[l|. In this case, the scale function reads L(t) = 1 + ut, 
where the constant u measures how fast the width of 
the wave packet changes (Berry and Balasz considered 
only the special case with u = 0). If u = 0, which is 
possible by preparing the initial wave function such 
that /3(0) vanishes, the intensity profile of the wave 
packet remains invariant while it experiences constant 
acceleration. In this case, xdt) = v^t + I/2agt^, where 
ug is the initial velocity and acceleration matches the 
constant og. Observe that acceleration is not due to the 
external force. If then the acceleration changes 

in time and accelerating wave packet expands/shrinks 
ballistically depending on the sign of u. 

Consider next the constant harmonic potential, 
= Wg. In this case, no physically acceptable 
self accelerating solution exists since x becomes singular 
when the scale function Lit) = cos(a;gt) equals to zero. 
If the harmonic potential is inverted, = — Wg [3^ . 

the self accelerating wave packet expands or contracts 
exponentially L(t) = depending on the initial 

condition /3(0) = 

Consider now time dependent harmonic potential. An 
interesting case is the non-accelerating wave packet 
solution in the presence of time dependent potential. 
Let us set ic = 0 in the coupled equations ( 00 ). We 
find that the acceleration of the wave packet is zero 
2^2 

if = —-2— where Vq = 9/2 og and Xq is 

9(a:g -I- Ugt)^ 

an arbitrary constant. This result contradicts with the 
Ehfenfest’s theorem, which states that the acceleration 
of the moving wavepacket is zero only for the free 
particle case. 

Instead of fixing the acceleration, one can also start 
with a pre-determined scale factor. As an example, 
one can study breathing wave packet by assuming 
Lit) = 1 -|- esin(r2t), where e is a small parameter and 
the constant LI is the angular frequency. We numeri¬ 
cally find that acceleration oscillates between positive 
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FIG. 1: The accelerations of wave packets in units of mmj^ 
when = —2(1 + 0.2cos(107rt)). We take ao = O.lmmjs^ 
and a;c(0) = 0.5/im. The thick curve is plotted according to 
the Ehrenfest’s theorem. The thin curve is the acceleration 
for the self-accelerating Airy wave. 


and negative values with growing consecutive peaks. 
Finally, let us study the evolution of wave packets 
for a given trap frequency. As an example, suppose 
oj^ changes periodically as = — Wq(1 -I- ecos(nt)). 
We solve the equations 00 with the initial con¬ 
ditions a;c(0) = 0.5/tm, i;c(0) = L(0) = 0. We take 
ao = O.lmm/s^, e = 0.2, = 2 Hz and = IOtt Hz. 

In the Fig-lU we compare the accelerations for the 
non-integrable self-accelerating wave packet and any 
integrable wave packet. The Ehrenfest’s theorem states 
that acceleration for integrable wave packets under 
such a time dependent trap oscillates non-periodically 
between negative and positive values with small ampli¬ 
tudes. However, the acceleration of our self-accelerating 
wave packet oscillates with growing peaks. 

Let us hnd the stationary solution of the equation (0. 
In the noninteracting limit. Go = 0, an exact analytic 
solution is available. It is given by 


/ / f\ \ ^0 \ 1 f\ 

Ipix ) = \ Ai{{ )3X ) , 

Ti 


( 7 ) 


where A is a constant and Ai{x) is the Airy function. 

The Airy function is not square integrable and decays 

exponentially for large positive x but oscillates for large 

t / 2.Tn^ ao , i 

negative x . One can see that the constant (—^—)3 

denotes the inverse width of the main lobe of the initial 
wave function. Hence we conclude that the constant oq 
is practically determined by the width of the main lobe 
of the Airy function. Therefore acceleration of the Airy 
wave can be changed by varying the initial width of the 
Airy function. In this paper, we show that self accel¬ 
erating waves have unique feature that the acceleration 
depends on both the external potential and the form of 
the wave packet. 

The Airy wave is a stationary solution in the accelerating 
frame. Transforming backward yields the accelerating so¬ 


I tcx, t=0) I ^ 



FIG. 2: The initial densities for self-accelerating waves when 
A = 5 X 10®. The thin and thick curves are for Go = 500 
(a typical order of magnitude for ultracold atoms in ID) 
and Go — 0, respectively. The densities are plotted with re¬ 
spect to the scaled coordinate X = {2m^aojh^)^^^x', where 
(2m^oo/?!.^)®^® = 10® m“®. 


lution in the lab. frame. 


tpix, t) = exp 



^c) + ^{x 



X 



2m? ao 
T? 



( 8 ) 


This is the exact analytical solution in the noninteract¬ 
ing limit. The wave function is not square integrable and 
does not obey the Ehrenfest’s theorem since the center of 
mass of the non-integrable Airy function cannot be de¬ 
fined. To this end, let us mention symmetry properties 
of self accelerating waves. The corresponding Hamilto¬ 
nian for our system remains invariant under both the 
parity V operation and time reversal T operation pro¬ 
vided that = =pw(t). However, the Airy function 

is not an even or odd function of position. Moreover, the 
time dependent function Xc(t) and the phase A(t) are not 
in general symmetric under time reversal. One can see 
from the solution (jS]) that V and T symmetries are spon¬ 
taneously broken for self-accelerating waves. 

So far we have considered the non-interacting limit 
where an exact analytical solution is available. Let us 
now study self-accelerating waves for the interacting case. 
It is worth saying that the relations (5,6) are not affected 
by the nonlinear interaction. Instead, the nonlinear inter¬ 
action changes the density profile of the self-accelerating 
wave. We now aim at finding the solution of 0. Since 
no exact solution is available, we solve it numerically. 
Let us first set the boundary conditions. The nonlin¬ 
ear interaction is dominant around the main lobe and 
gets weaker away from the main lobe. Since the density 
and consequently the nonlinear interaction go to zero as 
x' —>■ oo, the solution of 0 asymptotically coincides with 
the Airy function at large positive x . Therefore we set 
the boundary condition -ip ~ \Ai{(2m‘^ao/h'^y^^x') and 
dx'pJ Xdx'Ai{(2m'^ao/f?y^^x') at large x' to solve the 
second Painleve equation numerically. Before going fur- 
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ther, we remark that the constant A is a nontrivial de¬ 
gree of freedom in the interacting case. Increasing Gq at 
fixed ao and A increases the peak value of the main lobe 
slightly. If Go is very close to a critical value, Gc, then 
the peak value changes drastically. At exactly the critical 
2/2TO^ao, 


value, Go = Gc = A (- 




the peak value goes to 


infinity. Therefore the solution of the equation ([S|) is not 
bounded if Go exceeds the critical value, Gc- 
The Fig-2 plots the density of the self-accelerating so¬ 
lution at A = 5 X 10^ for Gq = 500m“^ (thin curve) 
and Go = 0 (thick curve). As expected, the effect of the 
nonlinearity can be mostly seen on the main lobe since 
the nonlinear interaction is weak for large \x'\. The peak 
value of the main lobe is increased by positive nonlin¬ 
ear interactions and the wave packet is extended with 
respect to noninteracting case. Furthermore, the main 
lope is shifted to the left when compared to the linear 
one. Note that no normalization is made on these solu¬ 
tions since they are not integrable. To this end, we would 
like to emphasize that time evolution of the nonlinear so¬ 
lution is exactly the same as the Airy wave packet (8). 
In other words, the solution given in the Fig.2 trans¬ 
lates while expanding/contracting in time according to 
(2), where Xc{t) and L{t) are the same for both nonin¬ 
teracting and interacting cases. 

We have shown that non-integrable wave packets lead 


to surprising physical results. However, realizing non- 
integrable wave function is impossible in practice. In¬ 
troducing an exponential aperture function is one pos¬ 
sible way for the physical realization of such waves 
^ 'k = where the decay factor e > 0 is a posi¬ 

tive parameter to guarantee the convergence of the func¬ 
tion on the negative branch Q. The Fourier transform 
of the truncated Airy wave packet is proportional to 
g-efe^gi(?i^/6m^ao)fc^ [^ . In photomcs Community, the 
Airy wave was experimentally obtained by passing a 
Gaussian beam throu^ a phase mask that adds a cu¬ 
bic phase modulation Since the wave packet is trun¬ 
cated, acceleration maintains for a finite distance only. 
The decay factor e determines how long the wave packet 
maintains its acceleration. For e << 1, accelerating be¬ 
havior of the ideal Airy wave packet was shown to be 
observed Q- 

To summarize, we have presented self-accelerating mat¬ 
ter waves for the harmonically trapped condensate. We 
have derived formulas for self-acceleration and shown 
that the initial form of the wave function and an initial 
phase contribute the self-acceleration. We have discussed 
that nonlinear interaction has nothing to do with the self 
acceleration but changes the density profile of the self- 
accelerating wave. It is worth studying self-accelerating 
waves in an optical lattice potential. 
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